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Abstract 

Spin systems are one of the most promising candidates for quantum computation. At the same 
time control of a system's quantum state during time evolution is one of the actual problems. It is 
usually considered that to hold well-known resonance condition in magnetic resonance is sufficient 
to control spin system. But because of nonlinearity of the system, obstructions of control of 
system's quantum state may emerge. 

In particular quantum dynamics of two 1 /2 spin-qubit system in the optical cavity is studied in 
this work. The problem under study is a generalization of paradigmatic model for Cavity Quantum 
Electrodynamics of James-Cummings model in case of interacting spins. In this work it is shown 
that dynamics is chaotic when taking into account center-of-mass motion of the qubit and recoil 
effect. And besides even in case of zero detuning chaotic dynamics emerges in the system. It is also 
shown in this work that because of the chaotic dynamics the system execute irreversible transition 
from pure quantum-mechanical state to mixed one. Irreversibility in its turn is an obstacle for 
controlling state of quantum-mechanical system. 
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I. INTRODUCTION 



Cavity quantum electrodynamics (CQED) is a rapidly developing field of physics study- 
ing the interaction of atoms with photons in the high-finesse cavities y|, |4| . Interest 
to such a systems basically is caused by two facts: One of them is the possibility of more 
deep understanding of quantum dynamics of open systems. Second argument is a possibility 
of practical application in the field of quantum computing In particular CQED experi- 
ments implement a situation so simple that their results are of great importance for better 
understanding of fundamental postulates of quantum theory [fj. They are thus appropriate 
for tests of bask; quantum properties: quantum superposition [?| , complementarily or en- 



tanglement 



101, 



11 



121 ] . In the context of quantum information processing, the atom 



and cavity are long-lived qubits, and there mutual interaction provides a controllable entan- 
glement mechanism an essential requirement for quantum computing In general 
dissipation processes must be taken into account when discussing problems of CQED. In 
particular there are two dissipative channels for systems the atom may spontaneously emit 
onto modes other then preferred cavity mode, and photons may pass through the cavity 
output coupling mirror. But modern experiments in CQED have achieved str ong atom-field 
coupling for the strength of the coupling exceeds both decay processes I3I, Q, Q- If so, 
then problem is reduced to the Jaynes-Cummings (JC) Hamiltonian, which models the in- 
teraction of a single mode of an optical cavity having resonant frequency with a two level 
atom comprised of a ground and exited states Q. 



1* 



One of the most promising candidates for quantum computation is spin systems 



19 



201 ] . In 21] was considered a two-spin-qubit system interacting with bath spins via 



Heisenberg XY interaction. The authors of indicated work could show that the problem is 
reduced to study JC two spin model. It has turned out that dynamics is non-Markovian. But 
in most general case atom- radiation field interaction should involve not only the internal 
atomic transitions and field states but also the center-of-mass motion of the atom and recoil 
effect. The study of such a case is the aim of this work. The subject of our interest is the 
following: it is well known that for quantum computing exact control of the spins system is 
necessary. That is why zero detuning is a matter of interest. In 22J has been shown that even 
taking into account of recoil effect and center-of-mass motion for zero detuning, dynamics 
is regular and chaos emerges only, when detuning is non-zero. But what will happen in case 



2 



of modified two spin JS model, it is not clear for the present. This work is devoted to the 
study of this problem. The first part of this work is devoted to quasi-classical consideration. 
In the second part we shall try to give kinetic consideration of the phenomena. 



II. QUANTUM NONLINEAR RESONANCE 



As was noted in the introduction we would like to consider more general model proposed 



in 



22j. It is not difficult to note that the Hamiltonian of the system of our interest 2l| 



takes the form when taking into account center-of-mass motion and recoil effect 221 ] : 

' +H s + Hsb + H B) (1) 



H 



2m 



where ^- is a kinetic energy of two spin qubit system placed in the resonator. It is supposed 
that qubit is composed of two spin 1/2 atoms 21J. The spin part of the Hamiltonian has 
the form: 

H s = u (SZ + S z 2 ) + n(SfSz + SiS+), (2) 

where H = 1, uq is Zeeman frequency of the spins being in the field inside the resonator, Q 
is a constant of dipole interaction between the spins in frequency units. The third term in 
(1) presents itself spin 1/2 atoms interaction with resonator field: 



H 



SB 



-g cos(k f x)((S+ + S+)b + (£r + S^)fr 



(3) 



here go is amplitude value of the qubit-field coupling that depends on the position of qubit 
x inside a cavity. The last term in (1) is the Hamiltonian of the field: 



H B = Ufb + b, 

where Uf is the selected frequency of radiation field, kf is the wave number. 
Taking into account commutation relation between operators 231 ] : 



(4) 



[b,b + ] = l, [S Z J ± } = ±S ± , [s + s- 



2S K 



it is possible to obtain the following Heisenberg equation of motions: 

dx p 
dt m' 



dt 



-g k f sm(k f x)((S^b + + Sfb) + (Szb + + S+b)), 



dS + 

— - = iu St - UnSlSt + 2igS{b + coa(k f x), 

(JjL 

dS{ 



dt 



-iuoS 1 + 2iVtS{S 2 — 2igS{b cos(fc/ 



(1^ Z A A A A A A A A 

= -igcos{k f x){S^b + - 5+6) - zfi(5 + 5 2 ~ - Sf S 2 + ), 

r7S + a a a 

2 - iu S+ - 2inS z 2 St + 2^S'f& + cos(A; / x), 



dt 

= —iu S 2 + 2iQS 2 S'-f — 2igS 2 bcos(kfx), 



dS 2 



dt 

dS Z A a A A A A A A 

= -ig cos(k f x)(S 2 b + - S^6) - iQ(S^S^ - S 2 Sf), 
db + 



dt 
db 



iuifb — ig cos(kfx)(Sf + S 2 ), 



-iujft i + ig cos(kfx)(S 1 + S 2 



dt 

After going to the representation of interaction: 

b+(t) = e iuJ f% b{t) = e-^H, = e^S 



± 



and introducing new variables by means of quasy-classical averaging \22\ : 



7 a < P > , 1 f+ f , 1 f 

x = kf < x >, p = — , o x = - < tr + o >, o„ = — < o 

; ' y k f 2 ' y 2i 

k 2 — 

a=— 5 = ^- — , p = Q/g , r = g t. 

mg g 

Taking into account (6), (7) we obtain from (5): 



dx 
dr 



— = ap, 



^ = -2smx{{Sfb x + S\b y ) + (£%b x + S y 2 b v )), 
ar 

dQ x 

—1 = -5S y + 2S{by cos x - 2PS Z S V 2 , 
dr 

dSf 

dr 

dS{ 



- 2£j% cosx + 2/3SfSZ, 



2 cos x(Sfb x - S*b y ) + 2j3(S?S% - S^Sf), 



d^ x 

—L = SS y + 2S z b y cos x - 2f3S z S y , 
—1 = SS X - 2S x 2 b x cos x + 2(3S z 2 Sl, 



^ = 2 cos x(S y b x - S x b y ) + 2/3(515? - 5f5f), 
rf6 



— — = — cosx(5f + 5f ), 

^k = - C0SX (S x + S 2 c ). (8) 
It is readily seen that the equations (8) allows the following integrals of motion: 

Sh = (St 2 ) 2 + (ST, 2 ) 2 + (ST, 2 ) 2 , N = bl + bl + S? + 55, (9) 

2 

W = ^- + 2/3(5*5* + 5?5f) - 2 cosx^*^ + 5 X \) + (5f 6* + S y 2 b y )) - 5(S Z + 5 2 z ). 
Introducing the new variables: 

U, = 2(S*b x + Sfby), U 2 = 2(S x b x + S y b y ), 

v x = 2{b y S x - b x S y ), v 2 = 2{b y S x - b x S y 2 ), 

g = (S x S y - S y S x ), f = (S X S X + S y S y ). (10) 

Taking into account the new variables (10) and integrals of motion (9), the set of equation 
(8) can be rewritten in more compact form: 

dx 

— = ap, 
dr 

^ = -2smx(U 1 + U 2 ), 
dr 

dS? 

— — = —2ui cosx + 2pg, 
dr 



dr 
dU x 

dr 
dU 2 

dr 



dS z 2 

- — 2v 2 cosx — 2pg, 



5ui + 2(3Sfv 2 — 2gcosx, 
5is 2 + 2fiS 2 v\ + 2gcos:c, 



^1 = SU! + 2 cosx(5 2 - 3(S Z ) 2 + 2NS Z 1 - 2S Z S Z + f) - 2f3S z U 2 , 
dr 

-p. = -5U 2 + 2 cosx(5 2 - 3(5 2 z ) 2 + 2iV5 2 z - 2S Z S Z 2 + /) - 2f3S z U 1} 



^ = cosx(5^ 2 - SZUJ - 2f3S{(Sl - (S z 2 ) 2 ) + 2^(Sf - (S*) 2 ), 

^ = (Sfo + Sfa) cos x. (11) 
ar 

By direct checking one can be convinced, that because of complex structure of the set (11), 
even for zero detuning 5 = 0, it is impossible to obtain analytical solution. Thus, unlike 



the problem studied in [22j, in case of qubit, taking into account of recoil effect and center- 
of-mass motion leads to nonintegrability of the problem even for zero detuning. Because of 
nonlinearity of the set (11) we expect to obtain chaotic solutions. If so, the state of qubit 
will not be possible to be controlled. 

We have integrated the set of equation (11) for the realistic values of parameters from 



the point of view of experiment [13j, [15| 5 = 0, a = 0.01, (3 = 0.5, Sf = ££=§. The 
results of numerical integration are presented on Fig.1,2. 

As is seen from Fig.l, the dynamics of the system even for zero detuning 5 = has chaotic 
form. 

The other parameters of the system have also chaotic spectrum (see Fig. 3). 
In order to be more convinced of dynamics to be chaotic, we have calculated even fractal 
dimension of the system. 

In order to calculate fractal dimension of the system's phase space we use the Grassberger- 



Procaccia algorithm 



24 



25] . The idea of this algorithm is the following. Let us suppose, 
we obtain an ensemble of state vectors Xi, % = 1,2, ...TV by numerical solving of the set of 
equations, corresponding to successive steps of integration of differential equations. Choosing 
small parameter e we can use our result for evolution of the following sum: 

1 - 

^) = iL-^F3T)g^-N--il)' ( 12 ) 

where 9 is a step function 



, x < 

6{x) = { (13) 
1 x ^ 

According to Grassberger-Procaccia algorithm, if we know C(e), we can estimate strange 



attractor's fractal dimension with the help of the following formula 24. |25| 



The numerical results are represented on Fig. 4. 
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FIG. 1: The graph of dependence of the system coordinates on time x(r). The graph is plotted 
for the following parameters x(0) = 1.6, p(0) = 9.1, 5f(0) = 5| = 0.863, £7i(0) = 0.000081, 
U 2 (0) = 0.000082, i/i(0) = 0.000083, u 2 (0) = 0.000084, g(0) = 0.0000845, /(0) = 0.0000846. As is 
seen from the plot trajectory has the chaotic form. 

II. As we have shown in the first part of the work the dynamics of the system is chaotic 
for certain values of parameters even for zero detuning. When considering the state of the 
system with quantum-statistical methods we shall neglect kinetic energy of the system and 
operator x will be regarded as classical chaotic variable x(t) , presented itself stochastic 
process. Condition of using this kind of approximation is the following: Acting on the 
system classical force is 

|F| = At" = \ V - H SB\-9oK f 
So, classical momentum transferred to the atom is AP m Atg Kf. Then influence of the 
atomic motion on the energy levels can be neglected if 1 2m ; < \Hsb\- 

ao Kj(At)* 

lhis means — i < 1. 
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FIG. 2: Fourier image of correlation function G p (t) = (P(t+r)P(t)) ,G p (u) = J dTexp(iuT)G p (T). 
Finite width of correlation function confirms the existence of chaos. The graph is plotted for the 
same values of the parameters as for Fig.l. 



Let us write Schrodinger equation of the system in interaction representation: 

.d\ip(t) > 



dt 



(15) 



where 



V = n(S+Sz + S^S+) + io f b + b - g cos(k f x)((S+ + S+)b + (Sf + 5 2 ")^ 



(16) 



is an interaction operator. 

Assume that at zero time t = 0, the system's wave function represents itself direct product 
of wave functions of atom \ip a tom > and \ip field > field: 

\lj)(t = 0) >= |VWn > Afield > ■ 



Here 



\lpatom >= COOICJO > +C l|01 > +C 10 |10 > +C U |11 >, 



(17) 
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FIG. 3: Fourier image of correlation function of variable Sf. The numerical vales of the parameters 
are analogous of that of Fig.l. 

\ipfield >= J2W n \n >, (18) 

where \n, m > is qubit's wave function. 

Because of interaction (16) the following transition between states are possible: 

|0,0,n+ 1 >^ |0, l,n >, |0,0,n+l ><-> |l,0,n>, (19) 

|0, 1, n ><-► |1, l,n - 1 >, |1, 0, n |1, 1, n - 1 >, |1, 0, n |0, 1, n > . (20) 

The transition (19) correspond to the transitions between energy states with changing 
number of photons and the transitions (20) correspond to inter spin transitions. On the 
basis of equations (19), (20) we shall search for the solution of equation (15) in the following 
form: 

,0,71+1 

|0, 0, n + 1 > + Co,i,n|0, 1, n > + 

n n 

+ J2 C i,o,n\h 0,n>+J2 Ci,i,n-i|l, 1, n - 1 > . (21) 

n n 

9 




FIG. 4: The graph of dependence of C(e) on e plotted using Grassberger-Proccacia algorithm 
for the values of the parameters analogous of that of Fig.l. A solid line corresponds to least- 
squares approximation of the results of date processing. The estimated fractal dimension is equal 



to D 



ln(C( £2 ))-;n(C(ei)) 
lne2— lnei 



2.2,(ei » 0.12, e 2 ~ 0.41, C(ei) » 0.5 • 10- 6 ,C(e 2 ) ~ 7.5 • 10~ 6 ). The 



numerical data obtained verify that the dynamics of the system is chaotic. We shall make use of 
this fact in the second part of this work where quantum-statistical description will be used for the 
study of the systems dynamic without use of quasi-classical methods. 



Taking into account equations (15)-(21), we obtain the following equations for coefficients 
of resolution: 



dt 

■ dCnn 



- = ttC 0fi , n +l ~ 90 COs(k f x(t))y/n+ r(C lj0 , n + C ,l,n), 



- = nCi,i jn _i - g cos(k f x(t))^(Ci fi ,n + Co,i,n) 



dt 

i = HCi.o.n - 9o cos(k f x(t))(C li i tTl - 1 y/n + C , ,n+iVw + 1), 

.dCiQ, 
dt 



^Co,i,n - go cos(kfx(t))(C lilin - iy /n + C 0fi , n +iVn + 1). 



(22) 
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In the set of equations (22) let us pass to the new variables: 



A{t) = C 1>0 , n + Co,i,n, B[t) = Vn + lC , ,n+i + VnCi,i,n-i. (23) 
Taking into account (23), the set (22) takes the following form: 

M{t) 



dt 

dB(t) 



QA(t) - 2g cos(k f x(t))B(t), 



QB(t) - (2n + l)g cos(k f x(t))A(t). (24) 



dt 

If we introduce the new notations: 

B'(t) = V2B(t), A 1 {t) = y/2n + lA(t), u(t) = V2~V2n + lg cos(k f x(t)) (25) 
and after that: 

C(t)=A'(t) + B'(t), D(t)=A'(t)-B'(t). (26) 

Taking into account equations (25) and (26), the set of equations (24) takes the simpler 
form: 

= QC{t) -u{t)c{t), 

(JjV 

It is readily seen that the solutions of the set (27) have the following form: 

C(t) = Jo ; D ^ = C2e int e -if^(t')dt'_ ^ 

Let us introduce the notations for the functionals: 

Q[w(t)]=e i fi u W* f , (29) 

Q*[u>(t)] = Q^Ht)] = e-ifXW, (30) 

Taking into account (27)- (30), the solutions of the set (24) takes the form: 

A(t) = -ff2= = Cl e- iQt Q[u(t)} + ° 2 e-^Q-^jt)], (31) 
V2n + 1 2V2n + 1 2y/2n + 1 

= = ^e-^QK*)] - ^ e -«*g-i[ w (t)]. (32) 
Taking into account (31), (32) and (23) we obtain: 

Ci,o,„ + C 0l i, n = - 3— v e- mt QHt)] + Sl— e-^Q-^t)], (33) 
2\/2n + 1 2\/2n + 1 

11 



x/^TTCo,o,„+i + v^Ci,i,n-i = -^e^g^(t)] - -9?-e- iUt Q- l [uj{t)\. (34) 

The equations (33) and (34) are the conditions to determine time dependence of the coef- 
ficients of the functions (21). But for determination of four coefficients we need two more 
conditions. The third condition for determination of coefficients C\fl^ n (t) and Co,i in (t) is 
easily obtained from equation (22) and has the following form: 

■ d(C ,l,n-g 1 ,0,n) = _ n{Co i n _ Ci Q ny (35) 

From this we have: 

Co,i,n — Ci,o,n = Cse lQt . (36) 
In order to obtain the last fourth condition, we introduce the notation: 



Vn+TC , ,n+i - v^Ci,i,n-i = F(t). (37) 



Then taking into account (22) we obtain for F(t): 

.dF{t) 



dt 

The solution (38) has the form: 



QF(t) - g cos(k f x(t))A(t). (38) 



t 

iQt 



■n -%\it /• 



iC 2 e~ 



iQt 



o 
t 



J Lo(t')Q^[io(t')}dt' + C 4 e~ lUt . (39) 



2^2{2n + 1 

u 

For further simplification of equation (39) consider the expression: 



i fu(t")dt" 

u(t')Q[u{t')]dt' = I u(t')e ° dt' (40) 





and let us introduce the notation: 

t' 

n (O = J oj{t")dt". (4i) 



Then it is readily seen that: 

t t 
u{t')Q[uj{t')]dt' = [ dQo(t')e iUoit,) = -t(e inoit) - 1). (42) 



o o 
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By analogy with previous one: 

u{t')Q^[uo{t')}dt' = I dn (t')e- inoit,) = i(e~ in ° {t) - 1). (43) 



Taking into account (42), (43), the expression (39) takes the form: 

f~i p —iVlt (~i —iflt 

m = Vl^TI) WMi)l " l) ~ vfe) (e " 1Kf)1 - 11 + C '- e " m (44) 

Taking into consideration (33), (34), (36) and (44) we can yet write down the set of four 
algebraic equations for the coefficients of wave function (21): 

2y2n + 1 2y2n + 1 

Co,i,n — Ci,o,n = Cse tQt , 



sjn + lC 0) o,n+l - V^CV.n-l = 

(Q[co(t)} - 1) - —J- -(Q-'Mt)} - 1) + C 4 e~ iQt . (45) 



2^(2^+1) 2^(2^ + 1 

Here the coefficients of integration are connected with the initial conditions via the rela- 
tions: 

d = v / 2^TT(C ,i,n(0) + C 1An (0)) + V2(V^+TC Ofi>n+1 {0) + v^7i, llft -i(0)), 



C 2 = v / 2^ + T(Co,i,„(0) + C lj0 ,„(0)) - >/2(Vn+TCo >0 ,» +1 (0) + >AIC7i,i, n _i(0)), 

C3 = Co,l,n — Cl,0,n) 

C4 = C ,o,n+i\/w + 1 - Ci Xn -i^/n. (46) 

by solving the set of equations (45), it is possible to determine time dependence of wave 
function (21) and by means of this to determine quantum state of qubit: 

^ +l(t) = 2^(2. + 1) " 2^(2. + 1) 

I ( ^ Cl + Ca ) e -<nt 

2v^TI 4 v / 2 v / ^T^(2n + 1) 4 v / 2v / ^TI(2ra + 1) 



C ie - int Jn ^ ... C 2 e~ int Jn „ u ... 
13 



^4 Cl C*2 n 

)e , 



2 v / ^Tl 4 v / 2v / ^TT(2n + 1) 4 v / 2v / nTT(2n + 1^ 

4\/2n + 1 4v2ra + 1 2 

Co,i,»W = , fef ^)] + , Af ^Wl " (47) 
4v2n + 1 4v2n +1 2 

As is seen from (47), time dependence of the coefficients of wave function (21) describing 
quantum state of qubit is determined by the functional: 

Q[u(t)]=e i fi u V ) * f , (48) 

where 

uj{t) = \/2{2n+ l)g cos(k f x(t)). (49) 

As is seen from (49) time dependence of quantum state depends on x(t). Thus in order 
to determine qubit's state, it is necessary to know the coordinate of the, system as explicit 
function of time x(t). But on the other hand as we have showed in the first part of the work, 
because of the dynamics to be chaotic x(t) may be considered as classical chaotic process. 
In this case for determination of the system's state it is necessary to average the functional 
(48) by all realizations of stochastic variable x(t). For this, we represent stochastic average 
of functional (48) in the form of the following continual integral: 

< Q[u(t)] >=< exp{i / Lo(t')dt') >— lim / duo N . . . du>xexp(i ^ u k At k )P N (u), (50) 

where 

N 1 

P N (u) = (2ir)- N / d\x . . . d\ N exp[-i ^ X k uj k }exp[-- C kkl \ k \ k >] (51) 

J k=i kk' 

is Fourier image of distribution function, At k = t (fc) - t (0) = 0, t (7V) = t. 

It is readily seen that by taking into account (51), the expression (50) can be rewritten 
in the following form: 

JV 



/ du)N ■ ■ ■ duiexp(i u k At k )P^ 
J k=i 



1 N 1 f 

d\i . . . d\ N exp{--^2c kk/ \ k \ kl ]Y\_7^ / exp[iuj k (At k - A fc )] 

kk' k=l 
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= I d\ x . . . d\ k 5(\ x - Atx)5(A 2 - At 2 ) . . . 6(X N - At N )exp[~ ^ C kk i\ k \ k i}. (52) 

kk' 

Taking into account (52) for statistically averaged functional we obtain: 
< Q[uj(t)} >= lim exp[-l-y C{t {k) ^ k>) )At k At k '} = 



kk' 
t t 



= exp(~ J dt' J dt"C(t\t")). (53) 
o o 

For random processes C(t',t") = C(t' — t"). Then introducing the new variables: t' — t" = t, 
t' + t" = £, and assuming that correlation function has Gaussian form 
C(r) =< uo{t + t)uj{t) >= e~ a ° T2 , finally from (53) we obtain: 

< Q[u>(t)] >= expi—.f^-Erfit^)}, (54) 
2 V a 



where Erf{. . .) is error function 26] . 



Assume that at zero time the system was in the state: 

1^(0) >= \ljjatom > ®\*Pfield >, (55) 

where 

l^fieid >= ^W n \n > . (56) 

n 

Comparing (55), (56) with: 

>=^C ,o,n+i(0)|0,0,n+l >+^C r o I i 1 „(0)|0,0,n> + 

n n 

+ Yl CWWIl, 0, n > + J2 Ci,i,„-i(0)|l, 1, n - 1 > (57) 

n n 

it is possible to obtain the following relations for the initial conditions: 

CooW n+ i = CoOn+l(0), CoiW n = 6*01^(0), 

C 10 W n = CW(0), C 00 W n+1 = Coon + i(0). (58) 

Let us determine the values measured on experiment that are connected with population 
difference of levels: 

hi,oi = W(t,\ll>)-W(t,\01>), 
15 



/n,io = W(t, \11 >)-W(t,\10>), 
ho,oo = W(t,\10 >)-W(t,\00>), 
Toi tOO = W(t,\01 >)-W(t,\00>), 
hi,oo = W(t, |11 >) - W(t, |00 >) = -(/ 11>0 i + /n,io) + 7^01,00 + Voo), (59) 

where 

oo 

W(t,|ll>) = ^|C 1 , 1 , n „ 1 (t)| 2 , 

n=0 

oo 



W(t,\01>) = Y,\C0Xn(t)\ 2 , 

n=0 

oo 

W(t,\10 >) = ^|C 1A „(t)| 2 , 



n=0 



W(t,\00>) = J2\Co,o,n + i(t)\ 2 . (60) 



n=0 



For illustration let us calculate for example W(t, |1, >). Taking into account (47) and (54) 
we obtain: 

-.00 ^ oo ^ 

< W(t, |10 >) >= - V(C 10 W„ + C m W n ) 2 + - V — — (v / ^TTC , ,oW / n + i+ 
8 ^— ' 4 ^— ' in + 1 

n=0 n=0 



+ V^C^Wn^f + ^ ^(C li0 W n - C ,iWg 2 + < |10 >) > mt , (61) 

n=0 

where < W(t, |1,0 >) >j n t denote interference terms whose explicit forms are not brought 
here for the sake of brevity. The point is that interference terms contain terms of the 
following form: 

< Q- 2 [uj{t)} >, < Q 2 [uj(t)] >, < e^'Q^K*)] > . (62) 

These quantities, as well as (54), fall down quickly after the lapse of time. For example: 

< e 2im Q- l [uj{t)} >= e 2im exp{-- x [^-Erf{t v ^)). (63) 

2 V «o 

As is seen from (63), for time interval that is more then the time of correlation function of 
the random quantity uj{t) (49), t > 



C(r) =< u{t + t>(t) >= e~ aoT (64) 
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interferentional terms can be neglected in (61). Situation is analogous for other quantities 
as well from (59), (60). Thus we were able to prove that because of dynamics to be chaotic 
zeroing of interferentional terms occurs. This fact of zeroing of inerferentional terms has 
deep physical sense. This means that the system execute transition from pure quantum- 
mechanical state to mixed one [23j. Such a transition is irreversible, as information about 



the phase of the system is lost. Transition from 



the manifestations of quantum chaos 



28 



29 



Dure quantum state to mixed one is one of 



30, 



31 



321 ] . Formulae analogous to (61) can 



be obtained for other quantities (60) as well: 



oo 1 1 oo 1 

< W(t, |0, >) >= J2 4 f 2n + + 1) {CloWn + ° mWn)2 + 2 £ J^V^^+^ C ^ n+ i+ 

n=0 ^ ' n=0 ^ > 

oo 

+V^c lil w n _ 1 ) 2 + V n (v^Co,otVn + i - V^TTc^w^) 2 

oo ^ oo 

< W(t, |11 >) >= £ 4(2n n + 1) iCloWn + CmWn)2 + 9 E 7^ ! ZlV2(^+ Tc o,oW n+1 + 



n=0 



n=0 



(2n + 1) 



+V^C 1)1 W n _ 1 ) 2 + ]T 



n=0 



n+ 1 
(2n + 1) 



r(\/nC 0) oWn + i - v/nTTCiiW n _i) 



1 oo 

< 1 10 >) >=< |0, 1 >) >= - J2(C w W n + C 01 W n ) 2 + 

n=0 

7 ^^(vVTTCo.o^n+i + v^d^n-i) 2 + - y^(C 10 W n - C m W n f 
4 In + 1 4 z — ' 

n=0 



(65) 



where C o, C i, Cio, Cn quantities are populations of corresponding levels, W n describes 
state of the field. It is usually assumed that W n satisfy Gaussian distribution ~\: 



(n — n) 2 
An 2 



-.expl 



n- 



(66) 



As we noted above transition from pure state to mixed one is irreversible. In order this fact 
to be confirmed, let us calculate change of the system's entropy. 

Let us assume, that the system at zero time was in state Coo, n +i ■ I n this case the system's 



entropy according to 



33|], |34j is: 

4 

S{t = 0) = -J2 Pi \n Pi = 0, 



8=1 



as only one of the elements of density matrix p is nonzero: 

Pl {t = 0) = |C 00 , n+1 (0)| 2 = 1, p 2 (t = 0) = p 3 (t = 0) = p 4 {t = 0) = 0. 



(67) 



(68) 
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After the lapse of time that is more than the time of transition between the levels to ~ 
l/go,t > to the system has time to execute transition between levels. That is why probability 
to find system in other states will be nonzero: 

Cn^O, Cqi^O, C 10 ^0 t>t . (69) 

Despite of this fact to talk about probability of population of different states is early yet. 
The point is that in time interval: 

t <t<J— (70) 
V "o 

interferentional terms in equations (61), (65) are nonzero. Therefore the state of the system 
will be pure one. But unlike of the initial state (68), which is simple state, the state of the 
system in time interval (70) is superposition one. 

Superposition state is pure quantum mechanical state and only after zeroing of interfer- 
entional terms in (61) and (65) superposition state passes to mixed one. Such a transition 
occurs in times: 

V «o 

But in time interval (70) while the system is in pure superposition state, from the symmetry 
point of view, it is clear that the coefficient values(69) have to satisfy the following relation: 




t/^J ~ c n lt <t< 



~C 01 [t <t<J^)~C 10 (t <t<J^)~C. 



Taking into account normalization condition: 



l 



< W(t, \ij >) >= 1 (73) 

i,j=0 

and (72), from (65) we obtain: 

oo 

C 2 {^ W n + W n+l + W n-l)) = L ( 74 ) 
n=0 

Then taking into account the relation: 

<W(t,\0,l>)>=<W(t,\l,0>)> (75) 
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it is easy to obtain the condition from (65): 



< W(t, |11 >) > + < W(t, |00 >) >=< W(t, |01 >) > + < W(t, |10 >) > + 

oo 

(+ Y, so n ^ (V^C 00 W n+1 - V^TTCnW^) 2 . (76) 

The condition (76) means in its turn that at times (71) mixed states are formed in the 
system in which the levels: 

Pl =< W{t > v^7a|H >) >= a, P2 =< W(t > v/ttHOO >) >= b (77) 

are populated with more probability than the levels: 

p 3 =< W(t > v^TalOl >) >= p 4 =< W{t > y^/allO >) >= c, (78) 

where quantities a, b, c satisfy normalization condition: 

Tr(p) = a + b + 2c = l, a + b > 2c (79) 

Taking into account (67), (77), (78) and (79) it is easy to see that during evolution of the 
system from pure quantum-mechanical state (68) to mixed one (77), (78) increase of entropy 
occurs. 

AS = S(t> v/^7tt)-5(t = 0) = -(alna + 61n6 + 2clnc) > 0, 

< a < 1, < b < 1, 0<c<l (80) 

III. CONCLUSION 

Let sum up and analyze the results obtained in conclusion. 

The aim of this work was to study two 1/2 spin qubit system being subject to resonator 
field. Interest to such a systems is caused by the fact that they are the most perspective to 
be used in quantum computer. The question that came up is the following: by how much 
will be state of the system controllable and dynamics reversible? We have considered the 
most general case, when interaction of the system with field depends on coordinate of the 
system inside resonator. 

Contrary to generally accepted opinion, it has turned out that the absence of detuning 
between resonator field and frequency of the system does not guarantee reversibility of the 

19 



system's state. During evolution in time the system executes irreversible transition from pure 



quantum-mechanical state to mixed one. At the same time, the time needed for formation 



of mixed state t > y ir/ao is determined completely by the system's parameters a = — -. 



38 



One more peculiarity of the problem studied is the following. It is well known 35 



36, 



39 



37|, 



40| that for integrable quantum systems complete and fractional quantum revivals 



are typical 



35 



36 



37 



38 



39 



40j. In our case because of dynamics to be chaotic phase 



incursion occurs. This results in zeroing of interferentional terms and irreversible losing of 
information about the system's state. This guaranties the absence of quantum revivals for 
our system, the noted fact may be especially interesting for experimental investigation of 
the system under consideration. 
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